Anomalous Scaling from Controlled Closure in a Shell Model of Turbulence 
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We present a model of hydrodynamic turbulence for which the program of computing the scaling 
exponents from first principles can be developed in a controlled fashion. The model consists of N 
suitably coupled copies of the "Sabra" shell model of turbulence. The couplings are chosen to include 
two components: random and deterministic, with a relative importance that is characterized by a 
parameter called e. It is demonstrated, using numerical simulations of up to 25 copies and 28 shells 
that in the N — > oo limit but for < e < 1 this model exhibits correlation functions whose scaling 
exponents are anomalous. The theoretical calculation of the scaling exponents follows verbatim the 
closure procedure suggested recently for the Navier-Stokes problem, with the additional advantage 
that in the TV — > oo limit the parameter e can be used to regularize the closure procedure. The 
main result of this paper is a finite and closed set of scale-invariant equations for the 2nd and 3rd 
order statistical objects of the theory. This set of equations takes into account terms up to order e 4 
and neglects terms of order e 6 . Preliminary analysis of this set of equations indicates a K41 normal 
scaling at e = 0, with a birth of anomalous exponents at larger values of e, in agreement with the 
numerical simulations. 



It is almost impossible to write a paper in the area of fluid turbulence 
without drawing heavily on one or more papers of Bob Kraichnan, 
who usually introduced the issue and wrote one or more fundamental 
papers towards its resolution. The present contribution is not an 
exception. It uses freely the model that Bob has introduced in j^j. 
Bob, we are happy to dedicate this paper to you at the occasion of 
your 10th birthday. The interaction, intellectual stimulation and 
friendship that you gave us over the last decade are unforgettable, 
and we hope and trust that they will continue for years to come. 



I. INTRODUCTION 



In this paper we present a calculation scheme aimed at 
evaluating the scaling exponents that characterize corre- 
lation functions of turbulent fields. For the sake of clarity 
we consider shell models of turbulence rather than ana- 
lyze the Navier-Stokes equations. However we stress from 
the start that in principle all the steps provided here can 
be repeated in the context of Navier-Stokes turbulence. 

Shell models as well as Navier-Stokes turbulence pose 
an infinite hierarchy of dynamical equations for the re- 
order correlation functions. This hierarchy is linear in the 
correlation functions, and in the limit of infinite Reynolds 
number is also homogeneous. It was recently discov- 
ered that this hierarchy obeys a rescaling symmetry 
which stems from the rescaling symmetry of the Euler 
equation Q. This rescaling symmetry foliates the space 
of solutions into slices of different scaling exponents h of 
the velocity; these are referred to as /i-slices. On each 
/i-slice one finds "normal scaling" with the given value of 
h. The full solution is a linear combination of all the so- 
lutions on the /i-slices with non-universal weights which 
are determined by the forcing on the integral scale of 
turbulence. Different orders of the correlation functions 
are dominated by different /i-slices, and accordingly the 
full solution has anomalous scaling. The anomalous ex- 



ponents are expected to be universal. 

In trying to evaluate the scaling exponents appearing 
in this way from first principles, we proposed |^.^,^[ to 
truncate the hierarchy of equations, preserving the funda- 
mental rescaling symmetry that gives rise to anomalous 
scaling. Truncation is problematic; in turbulence there 
is no natural small parameter, and therefore any closure 
of an infinite hierarchy is uncontrolled. It is therefore 
worthwhile to introduce a 1-parameter family of mod- 
els, characterized by a parameter e e [0, 1], which shows 
normal scaling when e = and recovers the anomalous 
scaling of the original model when e = 1. 

Following an idea in Q , we will construct such a family, 
and show that the transition from normal to anomalous 
behavior occurs at some finite value of e > 0. We will 
use e as a small parameter to regularize the closure proce- 
dure; we will show that our closed equations are valid to 
0(e 4 ), whereas the neglected terms are of 0(e 6 ). We can 
improve the closure scheme systematically by including 
terms of 0(e 6 ), neglecting terms of 0(e 8 ), etc. 

A way to achieve this small parameter is to couple N 
copies of the same turbulent system, be it a shell model or 
the Navier-Stokes equations, and choose the coupling to 
have both a deterministic and a random contributions, 
with relative amplitudes e and \/l — e 2 . For e = 1 we 
loose the coupling between the copies, and recover the 
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initial anomalous problem for any value of N. For e = 
we will show that that in the limit N — > oo we get nor- 
mal scaling. Thus for some value of e > and for large 
enough N we expect to see the birth of anomalous scal- 
ing, hopefully in a perturbative fashion. The existence of 
this transition is the main discovery of this paper, and we 
study it analytically using the e-controlled closure proce- 
dure, and by direct simulations of the iV-copied model. 

In Section 2 we review briefly the "Sabra" shell model 
of turbulence 0] , and introduce the copy space with the 
appropriate coupling. We study the resulting model nu- 
merically in Sect. 3. In Section 4 we start the construc- 
tion of the theory. First we present the statistical objects 
and derive the hierarchy of evolution equations that they 
satisfy, exposing in particular their scale invariance. In 
Section 5 we discuss the closure which preserves the scale- 
invariance, demonstrating how the smallncss of e is used 
to control consecutive closure steps. Sections 6 and 7 
are devoted to analysis of the resulting e-controlled clo- 
sure, and the results are compared in Sect. 8 with the 
numerical findings of Sect. 3. The paper is on the whole 
rather technical, and for the sake of the casual reader we 
summarize in some detail the points of principle in the 
conclusion section which is numbered 9. 

II. THE (N, e)-SABRA MODEL WITH RANDOM 
COUPLINGS 



where the phases 9 n are restricted by the set of equations 
6 n -x + n = #n+l ■ (5) 

Choosing 6\ and 82 arbitrarily, 8 n is determined for all 
n > 2. Evidently, the physical results of the model must 
be independent of the choice of the phases 8\ and 82. 
In particular, the only non-zero correlation functions are 
those which are independent of the phases 9 n . The non- 
vanishing 2nd and 3rd order correlators are 

S 2 (k n ) = (Kl 2 ), (6) 
■S 3 (fc„) = Im(u„_iM„)j* +1 ) . (7) 

Note that Eq. (Q) respects additional phase symmetry 

u n -> -u* . (8) 

One of the consequences of this symmetry is that 

Re(u„_iu„u* +1 ) = , (9) 

explaining why we only need to consider the imaginary 
part in Eq. ^) . The symmetries of this model were se- 
lected explicitly in || to give rise to a small number of 
nonzero correlation functions, with the aim of simplify- 
ing the calculations presented in the later sections of this 
paper. 



A. The Sabra model 

The starting point is the Sabra shell model as intro- 
duced in H: 



du n (t) 
dt 



ak. 



■n+ 1 



U n+ iU n+ 2 + bk n U*-1«7J+1 



-ck n -i U n - 2 U n -i 



vkl u n + f n (t) . (1) 



Here u n refers to the amplitude associated with "wave- 
vector" k n , where the spacing in this reduced model is 
determined by k n = ko\ n ; A is the spacing parameter, v 
the "viscosity", f n (t) a random Gaussian force which is 
operating on the lowest shells. The parameters a, b and 
c are restricted by the requirement 

a + b + c = 0, (2) 

which guarantees the conservation of the "energy" 



n 



£ = 5>„(i)| 2 , 



(3) 



in the inviscid, unforced limit. 

The equations of the Sabra model are invariant under 
the phase transformation 



u n — > u n exp(i6»„) 



(4) 



B. (N, e)-generalization of the Sabra model 

The standard available procedures to generalize dy- 
namical systems to N copies involve real variables . In 
shell models in general, and in the Sabra model in par- 
ticular, the amplitudes u n are complex. Therefore we 
rewrite Eq. ([!]), following B, in terms of the real and 
imaginary parts u' n = He{u n } and u" = Im{u n }. Doing 
this we guarantee, after the generalization to N copies, 
that the N — > 1 limit coincides with the original model. 
The equations are: 



du' n (t) 
dt 



7a,n+l( U n+l U n+2 + U n+l U n+2) 



+76,n(-«n-l u n+l + + 7c,n-l«- 



2 "71-I 



l n-2 a n 



71- 1) 



— vk 



2 u' 



(10) 



d<(t) 
dt 



7a,n+l«+l<+2 + K+lK+2) 



+Jb,n( U n-l u 'n+l + <-l<+l) + Ic.n- 1 (-<-2 u «- 1 



+<- 2 <-i) -<<+/;'(*) 



where 



7a, 7i — o,k n , 76,n — bk n , 7c, 7i — ck n 



(11) 



(12) 
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These equations can be written more compactly using 
the following matrices: 



A +1 = 
C +1 ^ 



1 
1 

-1 
1 



-1 

1 

1 

1 



(13) 



which also may be written as A a a , a „ and C^ l<7 ,, with 
a , a' , a" = ±1. The first subscript a' denotes line (+1 
corresponds to the upper line), second subscript a" de- 
notes column (+1 corresponds to the left one). Clearly 



Denoting u n , 
as 



a: 



a: 



c° 



(14) 



-i = u' n , w n ,+i = it" we rewrite Eqs. (|l^.|TT|) 



dt 



K>Ai 



+7b,nV>n-l,<T f Un+l,<T f ') +Ca>cj> 



-v)x n tin, a ~t~ fn.cr 



(15) 



As usual, we adopt the convention that repeated dummy 
indices (here a' and a") are summed upon. 

Next consider N copies of Eq. |[|). The copies are 
indexed by i,j or £, and these indices take on values 
-J,..., + J, 2 J + 1 = N, N odd. The ith copy is de- 
noted .o-. Let be the coupling between copies, 
which will be chosen later. Equations (|l^) for a collection 
of copies are written as 



dv m 
dt 



it 



'-+l U n+l,cr> U n+2,a" 



+7b,n«„-l,cr' u n+l,cr" J ' <r'«r"7c, n-1 u n-2,a' a n~l.a' 

-vkl u® a + f® . (16) 

To proceed we note that the index i is defined modulo N, 
and introduce a Fourier transform in the "copy" space: 



e=-j 



(17) 



Note that the index a is also defined modulo N = 2J+1. 
It is convenient to consider values a within "the first Bril- 
louin zone" , i.e from — J to J. We will refer to it as the 
a-momentum. Since Un] a is real, 



U n,a ~ \ U n,a) = U n,a 

In "a-Fourier space" Eqs. ( |l6| ) read 

du2. 



(18) 



dt 



2- = ^$ q ^.t[A 



0,1 



x {^V"ba,n+l^ +1)CT ,< +2iCT »+7b,n^_ 1)CT ,< +1)CT «] 
+ C,»V-l^ 2 ,,«Lu») - vk 2 n <, CT + . (19) 



where A Q p is the Kronecker symbol: A Q Q = 1 and 
= for a 7^ /3. Observe that we use Greek indices 
for denoting component in a-Fourier space, and Latin in- 
dices for copies in the copy space. As a consequence of 
the discrete translation symmetry of the copy index [i] 
Eqs. ( |l9| ) conserve a-momentum modulo N at the non- 
linear vertex, as one can see explicitly in the above equa- 
tion. The coupling amplitudes $ q '^<t m these equations 
are the Fourier transforms of the coupling amplitudes 
Our choice of these amplitudes will be presented 
in the next subsection. 



C. Choice of coupling 



We choose the coupling amplitudes according to 



Vn 1 1 



(20) 



where \p Q '^ are quenched random phases chosen with 
uniform distribution of the phase, independently dis- 
tributed with zero average. The meaning of quenched 
randomness in the context of a direct numerical simula- 
tion is that we run Eq. ( |l9| ) with a given random choice 
of $, obtain results averaged over the randomness of the 
forcing /, and then rerun with a fresh random choice of 
$. Only at the end we average the statistical objects over 
the runs. It is believed S, and checked numerically in 
this work, that when the number TV — > oo the statisti- 
cal functions are self averaging, and the last average is 
unnceded. 

The couplings satisfy the following symmetries: 

V|>a,/3.7 — \|/«i7>0 (\J) Q </ 3 >7)* _ \$)-<x,-0,-~< 

^<x,P,1 _ _ (21) 

The first of these conditions stems from the identity of 
copies, leading to the invariance of the equations of mo- 
tion to an interchange of copies in the nonlinear term. 
The second one is the reality conditions, the third im- 
poses energy conservation in the inviscid, unforced limit. 
The requirement 



(22) 



guarantees that for N = 1 we recapture the original 
model at any e. Note that for e = we have the so 
called "Random Coupling Model" proposed in the con- 
text of the Navier-Stokes statistics by Kraichnan in . It 
was understood that in the limit N — > oo the direct 
interaction approximation (DIA) becomes exact. After 
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proper elimination of the sweeping effect (in the frame- 
work on the Lagrangian-history D1A JTcfl or by means 
of the Belinicher-L'vov approach p!l|) the Kolmogorov 
1941 (K41) scaling appear as an exact solution of the 
Random Coupling Model. The same is true for the shell 
models (8|^] in which the sweeping effect is absent by 
construction. For e = 1 the coupling coefficients in the 
a-Fourier space ( |20| ) are index- independent. This cor- 
responds to uncoupling the equations ([l6]) in the copy 
space, because in this case = SijSij. Thus for 

e = 1 we recover the original Sabra model with anoma- 
lous scaling || . We see that with the choice of the cou- 
plings (|20| ) the model interpolates continuously between 
the "Random Coupling Model" for e = whose scaling 
is normal K41 (at N — ► oo) and the Sabra model with 
anomalous scaling for e = 1. A model of this type was 
proposed in the context of the Navier-Stokes statistics by 
Kraichnan in M. The consequences for the perturbative 
theory in the large N and small e limits were considered 
byEyink§. 
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FIG. 1. Structure functions \og 2 (S p (k n )) vs k n in log-log 
plot for p=2,3,4, e = 0.8 and N = 25. 



III. NUMERICAL INVESTIGATIONS 
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We are interested in the N — > oo limit of Eqs. (|l9|) as 
a function of e. As we discussed above, for e = we have 
the random coupling model with normal K41 scaling in 
the limit N — > oo, while for e — 1 the copies become 
uncoupled and with a proper definition of the structure 
functions all the scaling exponents ( n are the same as 
the scaling exponent of the original Sabra model. We 
thus expect that for intermediate values of e we may find 
values of £ n that interpolate between K41 and the Sabra 
value. 

This is actually what was found. We measured the 
scaling exponent (2 of the second order structure function 
S 2 (k„) = F 2 {k n -0), F 2 (k n ;t) being defined by Eq. @. 
The exponent £2 have been calculated either with a linear 
fit in the two decade inertial range, from n = 4 to n = 19 
(see Fig. [l]) or by using the fitting procedure on the 
all structure function introduced in [jfj. The results are 
the same although there is a difference in the estimate 
of the error bars. The error bars in Fig. || are the ones 
obtained for the linear fit. In Fig. || one can see the plot 
of the value of the anomalous corrections to Kolmogorov 
scaling, 8(2 — C2 — 2/3, as function of 1/N for e = 0.8 
together with the same curve for e = and for e = 1 for 
N ranging from 5 to 25. 

The equations of motion ( |l9| ) with a = l,b = c = — 0.5, 
were integrated with the slaved Adams-Bashforth algo- 
rithm, viscosity v — 4 x 10 -9 , a time-step At = 10~ 5 . 
The forcing was subjected on the first two shells, chosen 
random Gaussian with zero average and with variances 
such that (72/171 = 0.7 (in order to minimize the input of 
helicity H). Averages were taken for a time equal to 
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FIG. 2. 5( 2 = Ca ~ 2 / 3 vs 1 / N for 6 = °- 8 (diamonds) and 
e — (circles) for N from 5 to 25. 

250 eddy turnover times for the case N = 1. The averag- 
ing times were decreased when the number of copies in- 
creased, taking into account the faster convergence times 
in these cases. 

The random couplings were chosen with a zero-mean, 
uniform probability for the random phase at the begin- 
ning of each simulations. Rigorously one should have 
taken averages over different runs with different couplings 
but we checked that self-averaging is already valid for 
N = 5, at least for e = 0.8 (very small random compo- 
nent in the couplings) and within our numerical preci- 
sion. For e = self-averaging occurs for large numbers 
of copies. 

In Fig. ^ it is clear then, while in the case for e = 
the corrections to Kolmogorov scaling goes to zero, as 
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was already checked for the equivalent spherical model 
||, for e = 0.8 and for e = 1 the corrections converge to 
a finite value which increases with e. 



IV. THE EQUATIONS OF MOTION OF THE 
STATISTICAL OBJECTS 



A. Multi-time correlation functions: definitions and 
symmetries 

The nonzero 2nd and 3rd order multi-time correlation 
functions of the collective variables are defined as 

F^ a ,{k n ;t-t')A a , a , = « CT (i)<MO> , (23) 

(24) 

+A a+a , +N , a „] ee (<_ 1)ff ,(i)< CT ,(tO<+i>(*")) ■ 

The conservation of the a-momentum which was dis- 
cussed above causes the 2nd order correlation function 
to be diagonal in a, a' , and puts a constraint a + a' — 
a" modTV on the 3rd order correlation. In (^3|) we wrote 
explicitly that in stationary conditions F 2 depends on the 
time difference t — t' only; it is clear also that F3 depends 
on two time differences, say t — t' and t — t". There are 
a few nonzero 4th order correlation functions, we present 
here a natural general definition: 

^4;cti ,0-2,0-3 ,(T4 (^«i ? ? i j ^1 : ^2: ^3) ti) A ai +ct 2 ,a 3 +ai 

ee ki^uz^uzu^uz:^)) • ( 25 ) 

As we see the correlation functions F n defined by (p3|)- 
( p5| ) depend on the indices a and a. This dependence is 
determined by symmetry considerations. The original 
Sabra model ([!]) is invariant under the transformation 
(||) . This leads to the following symmetry of the N, e- 
Sabra model (p6|): 



(26) 



Consequently, all the correlation functions must be in- 
variant under this transformation. For F 2 it means 
F£ t '" t<r ,(k n ;t) = era' F^ a ,{k n \t) and thus era' = 1, i.e. 
diagonality in a, a': 



F£,a : <j'(kri',t) — A a ^i F2 a {k n ] t) 



(27) 



The nonzero components of the .F3 tensor are those for 
which the product a a' a" = 1, namely: 



rria,a a j-ia.a a rpa 7 a a 

r 3;l, 1,1 ' -^3;1,-1,-1 > -^3;— 1,1,- 



1 ' -^3;— 1,-1,1 



(28) 



The nonzero components of F4 (|25|) are those in which 
the product o\ a 2 03 04 = 1. The corresponding list is 
obvious. 



The symmetry of the Sabra model under the phase 
transformation (Q, ||) is more restrictive. For the (N,e)- 
Sabra model it leads to an invariance under the transfor- 
mation: 



l n, + l 



cost 
sin( 



— sinf 
cos 9, 



a 

n,-l 
a 

n, + l 



(29) 



when 8 n + 6 n+ \ = 9 n+2 . The consequence of this for F 2 
is rather simple: the only invariant combination of the 
2nd order correlation function is the trace 



a=±l 



F?Jk n ;t) 



(30) 



Actually we can make an even stronger statement: the 
action of the random force on the first two shells results 
in the randomization of the phases 9% and 9 2 . Since 
the other phases are determined, they are also random, 
and satisfy (cos 2 9 n ) = (sin 2 6>„) = 1/2. Therefore 

F 2yCT {k n ]t) — ^F^ykn'ft). 

Moreover, the function F 2 (k n ;t) is a-independent 



F? a (k n -t) = \F 2 (k n -t) 



(31) 



This is equivalent to the statement that the correlation 
functions are diagonal in the copy space: 



(u [e] u [e ' ] * 



oc At 



(32) 



To see this note that the deterministic part of coupling in 
Eq. ( [l6| ) is diagonal in the copy space, i.e. oc eA^-A^. 
The "inter-copy" part of coupling (which is oc yl — e 2 ) 
is random (being a Fourier sum of the random phases 
tf^'T). Therefore every realization of Ws implies differ- 
ent phase relationships between the variables u n \a with 
different values oil. Therefore the only second-order cor- 
relation functions that survive the ^-ensemble averaging 
are those that are diagonal in £, £'. 

Finally we conclude that among the group (|23|) of 
2nd order correlation function there exist only one ob- 
ject, F2(k n ; t), which is invariant under all the symmetry 
transformations of the model and is independent of the 
indexes a and a. In terms of (|23|) it may be written as: 



F 2 (k 



OL,<7 



(33) 



In a similar but much more involved manner one may 



find the unique invariant object among the group (24) of 
3rd order correlation functions: 



1 



a,/3,7 cr.a' a" 



+ A^ +7+A r]F 3 ^;> (fc„; t, t', t") . (34) 



The easiest way to derive this expression is to compute 
a time derivative of the invariant object F 2 (pq) as is 



done in Appendix A 1. The only objects which may ap- 



pear in the result are invariant combinations of the 3rd 



(i 



order correlation functions and this is exactly the combi- 
nation . Computing a time derivative of F 3 (|34|) we 
find the following invariant combination of the 4th order 
group (|f): 



;*2>*3)*4)*l) • (35) 



x p /3, a, 

4; 0"2,CT3,(T4 



There are few other invariant combinations of the 4th 
order group (^5|), but only the combination (|3^) appears 
in the dynamics. 

Note that the normalization constants iV _1 in the defi- 
nitions (p3|,|34]) and A~ 2 in ( ^5| ) are chosen such that F 2 , 
F 3 and F 4 are of 0(1) . Moreover note that although the 
triple correlators ( |24|) a re complex functions the invari- 
ant combination F 3 (|34j) is real as a consequence of the 



reality condition for u a and of the sums over a and f3. 
In addition, all these definitions of correlation functions 
coincide with the original correlations of the Sabra model 
when N = 1. 



B. Hierarchical equations 

We now present the hierarchy of evolution equations 
for correlation functions up to the equation for F3, which 
contains F4 on its RHS. We will close the hierarchy by 
expressing F4 in terms of lower order objects, showing 
that what we neglect is of higher order in e. 

With the definitions given above, in the inviscid un- 
forced limit (i.e. in the bulk of the inerti al ra nge) the 
evolution equations derived in Appendix A 1 take the 
form: 



0_ 

di 



F 2 (k ni t) = 7 a ,„+iF 3 (A; ri+ i;0, t,t) + -fb, n F 3 {k n ; t, 0, t) + 7 C)n _iF 3 (fc n _i; i, f, 0) 



(36) 



F 3 (k n , t 1 ,t 2 ,t 3 ) — J a ,nF4(k ni k n+ i, k n , k n+ i- 1 t 1 ,ti,t2,t 3 ) 
~t~7f>,n— iFi{k n — 2i k n , k n , k n +i 

; tl,t l7 t 2 ,t 3 ) + 7 Cj „_ 2 i<4(fcn-3, &n-2, k n , k n+ i ] 1 1 , t\ , t 2 , t 3 ) , (37) 



■7T—F 3 (k n , ti, t 2 , t 3 ) — Ja,n+lF4(k n ^i, fc„+l, k n + 2l fc n +lj ti, t 2 , t 2 , t 3 ) 
Ot 2 



+"fb,nF4(k n -i,k n ^ 1 , k n+ i,k n+ i;ti,t 2 ,t 2 , t 3 ) + 7 Ci „_iF4(fc„_i, k n ^ 2 , fc„_i, fe„+i; t\, t 2 , t 2 , t 3 ) 

■7T—F 3 (k n , ti, t 2 , t 3 ) = 7a,n+2^4(fcn-l, k n , k n + 2 , k n + 3 ] t\, t 2 , t 3 , t 3 ) 

ot 3 

~^~^fb,n+iFi{k n —i, k n: k n , k n + 2 ; ti,t 2 , t 3 ,t 3 ) + 7 C! „F 4 (fc„_i, k n , k n , fe„_i; t-y, t 2 , i 3 , t 3 ) , 



This set of equations exhibits terms on the RHS com- 
ing from non-linear contributions only. The linear dis- 
sipative terms are negligible in the limit v — > 0: there 
is nothing in the dissipative terms that compensates the 
vanishing of the viscosity. In addition, the terms coming 
from the forcing appear explicitly only in the equations 
for the correlations calculated in k\ and k 2 , the forced 
shells. They vanish in the inertial range. 

The structure of the hierarchical Eqs. (|3^)-(]37|) and of 
various future equations becomes transparent if we adopt 
a graphical representation. The symbols used in the 
graphical representation are those introduced in the con- 
text of renormalized perturbation theory |0,p|, see Fig. . 
It is important to stress however that the present ap- 
proach is not perturbative, and does not generate any 
infinite series of diagrams. Accordingly the present the- 
ory does not suffer from the usual perturbative problems 
of unknown convergence properties of infinite series. 

The hierarchical Eqs. (|36|)-([37|) are exhibited in 
Fig. |]a. The consistent closure procedure (see below) 
requires the introduction of response functions and their 
equations of motion. The equations for all the response 



functions up to 4-point objects are presented analytically 
in Appendix |A 3| and graphically in Fig. ^b,c. 

C. Scaling invariance and ft-slices 

The Sabra and (N, e)-Sabra equations ([[]) and (|l6|) in 
the unforced, inviscid limit are invariant under the fol- 
lowing 2-parameter rescaling group lZ(h, p): 

K(h,p)k„ = k„/p, K(h, p)t = p^H , 

K(h, p)u n = p h u n , (38) 

where 1Z is the rescaling operator. Correspondingly the 
hierarchical equations for the correlation functions and 
the Green's functions display invariance under the same 
rescaling group lZ(h, p) when the statistical objects trans- 
form as follows: 

TZ(h, p)F n = p nh F n , TZ(h, p)G nA = P {n - 1)h ~ 3 G nA , 
K(h,p)G h2 =p-( h +VG h2 ,... (39) 
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a. Correlation functions: 



T, - AAAAAA i T. 




b. Greens' functions: 





6,3 = AAA^) ,. 



c. Vertices: 

Y ABC 

< < < < 



4-n,n 3,1 2,2 1,3 



lr , ~ ir ir ir ir ir , 
D = D , D , D , D D 

S-n." 1 4,1 3,2 2,3 ' 1,4 ' 0,5 



FIG. 3. The diagrammatic nota- 
tion of the basic objects, Panel a: The 
correlation functions T n . Panel b: 
the bare Green's function (thin line) 
and the dressed Green's functions 
Gn,m- A circle connecting n wavy 
lines stands for an nth order corre- 
lation function, and a circle with n 
wavy lines and m straight lines stands 
for a Green's function with n veloci- 
ties and m forcing. Panel c : the 3rd 
order vertices 7, and A, B, C; irre- 
ducible contributions of the 4th and 
5th order vertices D"_ n n and _Ds_ nin 
denoted as empty squares and pen- 
tagons. Double line tails stands ei- 
ther for a straight or wavy tail. 



a. Equations for correlators , 



a 1 




b. Equations for Greens' functions Q nAh : c. Equations for Green's functions Q h2th ■ 



ST? 



AAA = 



+ AAA 5 (t 12 ) , 



d ' 

<:// 





a 1 

(777 




FIG. 4. The symbolic representation of the first evolution equations in the hierarchy for correlation and Green's functions. 
The graphical notations are presented in Fig. 3 



Explicitly, the application of 1Z(h, p) to any statistical 
object reads: 

T^(h, p)F n (k m j k m i , k m " , • • • ;t, t',t",---) (40) 

__ „ (k m k m > km" „l-h+ n l-h + l J--h,n \ 

V p p p J 

The simple fact that the hierarchical equations are lin- 
ear equations has a crucial consequence for the theory: 
by direct substitution one can check that the rescaling 
symmetry is wider then (^9|) and includes an additional 
rt-independent function of h - 2(h), which at present 
seems arbitrary: 

K(h,p)F n =p nh+z ^F n , 
TZ(h,p)G nA =p (n -^ h - 3+z(h ^G nA , 
n(h,p)G ia = p-^ +z ^G 1 a 1 ... (41) 

Note that in principle the rescaling group for the hier- 
archies of the correlation and Green's function could be 
defined using two different functions 2\(h) and 22(h). 
However, the existence of an inhomogeneous term pro- 
portional to the bare Green's function (first line in 
Fig.^b) results in the identify of 2(h) in the two hier- 
archies. This point is explained in detail in pj. 



It is tempting to offer a physical explanation of the 
identity of the of Z(h). In some way it must be con- 
nected with the scaling exponent of the probability distri- 
bution function of subset of velocity configurations with 
a given h. In other words, if the probability to find in 
a turbulent ensemble a solution {u n } of the unforced 
equations in the inviscid limit which is Pt ({«„}), then 
the probability to find a rescaled solution 7Z(h, p){u n } 
is p z ( h > Ph(1Z(h, p){u n }). It should be stressed however 
that the issue of the rescaling properties of the ergodic 
measure of our model is far from being exhausted by this 
comment, and a lot of further research is needed to so- 
lidify and interpret it further. We leave this interesting 
issue for future research. 

We reiterate that at this point p and h are free param- 
eters, and Z(h) is an n-independent arbitrary function of 
h. We will see later that this freedom is restricted by the 
equations: although p remains arbitrary, h belongs to a 
natural interval, and the freedom of 2(h) is removed by 
a solvability condition for the hierarchy of equations. 

As stressed before M, the existence of the rescaling 
group TZ(h, p) for the linear set of equations suggests the 
existence of particular solutions that are characterized 
by h. In other words, we consider solutions denoted as 
F n -h and G m . n -h that solve equations on an h-slice, for 
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example 

Q 

-7rrF2-h(k n ,ti,t2 
ati 



lb,nF^-h{k n ; ti, f2, fi) + 7 Ci „„iF 3; / l (A:„_i; ti, ti, i 2 ) • 



(42) 



where the rescaling property of a particular solution F n h 
is 



H{h, p)F Ujh 



_ nh+Z(h) 



Fn,h(krn i tx ; • • • j t<n) 



(43) 



The general solution F n of the hierarchical equations that 
we are interested in is naturally obtained as the weighted 
integral 



F n ({k m },{t})= / dn(h)F n , h ({k m },{t}), 



(44) 



where {k m } — k m i, k m », ■ ■ -and {t} — t', t" , ■ ■ ■ are the 
sets of k- and t- arguments of correlation functions F n . 
The weight dp{h)/dh may be, in principle, found from 
the boundary conditions, by matching the general solu- 
tion in the inertial interval with non-universal part of F n 
in the energy containing interval. 

The representation (|44|) follows rigorously from the 
rescaling symmetry of the problem and the linearity of 
the hierarchical equations. It offers a very useful connec- 
tion to the multi-fractal model for anomalous scaling in 
turbulence that was introduced phcnomcnologically for 
simultaneous fully fused objects in In [0 it was re- 
alized that the multifractal model follows naturally from 
the equations of motion for the fully unfused multi-time 
correlation functions. Introduce the following dimension- 
less objects: 



F n , h ({K m };{T}) = — 



nh+Z(h) 



F„, h ({k m };{t}), 



(45) 



where k m is the geometric mean of the {fc m }, vectors 
involved in the correlation function: k mi , k m2 , • • •. Di- 
mensionless k - and i-arguments k and r are defined as 



U 



k 



t(rrij, h) 



(46) 



where t(m, h) is a characteristic time on the m-shell for 
the ft,-slice is 



1 /k \i-h 
k a u V k ) 

Finally the integral (E3) is written as 



(47) 



F n ({k},{t}) = u% / dn(h) 



k 



nh+Z(h) _ 



F n , h {{n},{r}) . 

(48) 



It is easy to recognize that this formula is a natu- 
ral (multi-time, multi-point) generalization of the multi- 
fractal representation of structure functions if Z(h) — 
3 — D(h) where D(h) is the fractal dimension of the set 
of points that support an Holder exponent h Ej. The 
function Z(h) is related to the scaling exponents of the 
structure functions £ n via the saddle point requirement 



min \nh 

h 



Z(h)] 



(49) 



This identification stems from the fact that the integral 
in (^8|) is computed in the limit (ko/k m ) ~ * o° via the 
steepest decent method. Neglecting logarithmic correc- 
tions one find that F n oc (ko/k m p n . 



V. ^-COVARIANT CONTROLLED CLOSURE 

In [0 we introduced a closure scheme that exploits the 
rescaling symmetry of the hierarchy of equation on h- 
slices. This closure turns out to be somewhat different 
from traditional closures. Since the first attempts in the 
40 's it has been customary to close the set of equations 
in the most natural way, that is by expressing high order 
correlation functions in terms of lower order ones. Gen- 
erally such procedures do not preserve the fundamental 
rescaling symmetry (|38|). For example, if one tries to es- 
timate F4 as F$, it would require Z(h) — 0. We propose 
that attempts to close the hierarchy without respecting 
the rescaling symmetry (with Z(h) ^ 0) necessarily lead 
to fundamental mistakes of this type, ruining the possi- 
bility of finding anomalous exponents. 

The minimal requirements for an acceptable closure 
should be that realizability conditions are guaranteed: we 
cannot have negative PDFs and negative energy spectra. 
In addition, we want to preserve the Galilean and rescal- 
ing symmetries of the original equations. We therefore 
propose a closure scheme in which both these require- 
ments are fulfilled. Realizability is guaranteed because 
in our model we have small parameters, i.e. 1/N and e 
that are used to control the closure, making the neglected 
terms vanishingly small compared with those retained in 
the limit 1/N — > and e — * 0. The second requirement is 
met by expressing higher order correlation (fourth order 
quantities in the lowest order closure) in term of lower 
order ones in such a way that the rescaling symmetry 
( ^l|) in an /i-slice is preserved. In particular all our clo- 
sure steps are Z(h) covariant, in the sense that all the 
terms in the resulting equations have the same Z(h) de- 
pendence. As a consequence the equations are neutral to 
power counting, and Z(h) is computed from coefficients 
rather than powers. In fact Z(h) is found from solvability 
conditions as is demonstrated below. 

The procedure developed in the rest of this paper is 
a closure at the level of F^h- We are going to represent 
F^h in terms of the lower order quantities, the two point 
quantities F^-h and Gi.i-.h, and the three point quantities 
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F^,h, Gi,2-h and (?2,i;/i- There is a part of F^-h that can- 
ned be represented in terms of these lower order quanti- 
ties, but we will show that this part is of lower order in e, 
in fact of 0(e 6 ) whereas the retained terms are of 0(e 4 ). 
In order to show this explicitly we will have to consider 
the closure at the level of F 5 -h, representing F 5 -h in terms 
of all the lower order 2-point, 3-point and 4-point objects. 

In representing higher order quantities in terms of 
lower order ones there are many different combinations 
of low order quantities which appear in any given higher 
order one. In fact all the topologically possible combi- 
nations are allowed. To facilitate the analysis it is nat- 
ural to represent these by diagrams. In making three 
point objects from two point ones, there appear three 
different ways of combining the two point objects. These 
are represented by 'vertices" which are denoted below as 
Ah, Bh and Ch- In perturbation theories such vertices 
are commonly represented in terms of infinite series. In 
our approach (as in a fully resumed diagrammatic the- 
ory) these vertices are solved for exactly in terms of the 
three-point objects on an /i-slicc. It turns out, once the 
theory is written down, that these vertices appear always 
in the same combinations with the 2-point Green's func- 
tions Gii-h which never appears alone. We thus work in 
terms of three newly defined objects that we call "trip- 
lices". The A-triplex is the vertex Ah linked to Gx,i-,hi 
the i3-triplex is the vertex Bh linked to two Gi,i ; fc's, and 
the C-triplex is the vertex Ch linked to three Gi^i^'s. 
We reiterate that at no point do we need to sum infinite 
series. 

At the end of our procedure we will pose equations 
of motion for the F 2 -h and the three triplices, which are 
closed upon themselves and are valid to 0(e 4 ). A simpli- 
fied version which includes only one A-triplex is shown in 
Fig.|^. This simplified set of equations is analyzed in the 
sequel, and we show that as a function of e it predicts a 
transition from normal scaling to anomalous scaling. The 
reader who is not fascinated by the technical details can 
skip Sect. VI, and obser ve the appearance of solvability 
conditions in Sect. 



VIII 



VI. 



TECHNICAL DETAILS OF DERIVING THE 
CLOSURE EQUATIONS 



A. Definition of the triplex objects Ta-hi Ts;hi To ; h 

As a preparatory step for the closure procedure we in- 
troduce three types of triplices: A-triplex Ta-,Hi B-triplex 
Tb ji and C-triplex Tc-h- Triplices are functions of three 
k- and three t arguments, say k, k' , k" and t, t' , t". 
As any three point objects of the Sabra-shell model they 
differ from zero only for three consecutive shell numbers, 
say m— 1, m and m + 1. The C-triplex is invariant under 
all permutations of points, therefore it is enough to 



One triplex approximation : 



JL AfWW- 
at i 



+ — -i > vAAA 

2 %AAN 




FIG. 5. The hierarchy of equations on /i-slice after the 
level-4 closure. For the beginning we did not consider here 
contributions of B- and C- triplices. This allows to reduce a 
system of four equations to the system of just two equations 
for F2,h (two top lines) and for A-triplex (two bottom lines). 
Equation for yl-triplex follows for equation in Fig. 8 after dis- 
regarding the empty square term and multiplying the result 
on the right by two inverse Green's' functions, fxlosedl 

mention only one of the fc-arguments, corresponding, for 
example to the middle shell k m . So, we have Tc.h = 
Tc;h(k m ;t,t' ,t"). The convention here is that the first 
time argument t corresponds to the smallest shell m — 1, 
the second t' to the middle shell m and the last t" to 
the m + 1-shell. The same convention is chosen for the 
A- and i?-triplices, the only difference is that these two 
are invariant only under permutation of two points, the 
third one is special. Therefore we have three different A- 
and i?-triplices (and one C-triplex) denoted as: 



T 



(+i) T (-i) 



A-h 



A:h 



T 



(0) 



A:h ' 



T 



(+i) T (~i) 



B;h 



B:h 



T 



(0) 



B:h ' 



(50) 



The superscript (+1) means here that the /c-argument of 
the special shell is k m +i, (—1) corresponds to fc m _i and 
(0) - to k m . These objects are defined in terms of the 
following seven three point functions (with the same con- 
vention concerning the arguments, see their definitions 



dA4j)-(|A7D: 



^1,2:71 ' 



r (o) 



G 



(+i) r (-i) 

2,1-h i ^2,l\h ' 



G 2 ,i:h .) 



(51) 



Three A-triplices Ta-h are defined in terms of the three 
corresponding (with the same superscript) Green's func- 
tions G\fl;hi which describe the response of the velocity 
v rrll to two f orcin g factors on the shells 1112 and raz de- 
fined by Eq. (A7). In terms of G\ t 2;h we have: 
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poo pOO 

Gfl h (k m ,t 1 ,t 2 ,t 3 ) = dt[ d^T^(*v n ;ti,t2,OGi,i;fc(*v n -i;ti-ti)Gi,iifc(fem+i;t3-*3) • (52) 

Jtl Jt 3 

Similar equations serve to define the triplices T A } h and T A .^ simply by replacing on the LHS the appropriate response 
functions and by the required changes of the time integration on the RHS: 



(53) 



poo poo 

G { ^l{k m ,t u t 2 ,h) = / dt[ ^ri+ 1) (fc m ;i' 1 ,t' 2 ,t3)G 1 , 1;ft (fc„ l _ 1 ;t' 1 -t 1 )G 1 , 1;ft (fc m ;t' 2 -t 2 ), 

Jtx Jt 2 

/>oo />oo 

G l 2 . h (k m , ti, h, t 3 ) = / dt 2 / dt 3 T A . h (k m ; t\, t 2 , t 3 )Gi t \-h{k m ; t 2 - t 2 )G\^-h{k m+ i\ t 3 — t 3 ) . (54) 
Jt 2 Jt 3 



a. Green's functions <? 1>2 . A : 




c. 3'rd order correlator 




FIG. 6. Exact representation of the 
Green's functions G2,v,h and Gi,2;h 
(Panels a and b) and correlation func- 
tion F3-h in terms of the vertex A and 
the Green's' functions Gi,i. In our 
approach these serve as the definitions 
of triplices Ta-,h, T B -h, Tc-h, on an h 
slice via F 2 - h and G 2 ,i.h, G 1;2l h, F 3 . h . 



Now let us have a look in Fig. 6b at the representation of the Green's functions G 2 ,i in terms of the two triplices 
A, B, the correlation function F 2 and the Green's' functions G\ t \. This serves to define the B-triplcx. The explicit 
analytic definition of T B °) h triplex is: 

oo oo oo 

G 2 °l. h (km,ti,t 2 ,t 3 ) = J dt 2 j J dt' 3 T^ h (k m \ti,t 2 ,t' 3 )F 2 -h(k m+ i]t' 3 — t 3 )Gi t i-h(k m \t 2 — t 2 ) + J dt[ 

t 2 

oo 

x rit fc 1) (^;^i,^,t3)G M , i fc(*^;^-t 2 )F2;fc(«5m-i;ti-ti)} + J dt' 2 T B % (k n ■ h , t' 2 , t 3 )Gi,i; fc {km-! ; t' 2 - t 2 ) (55) 



Similar equations for T B ~^ and T B x h can be derived by using G 2 ~^ h and G 2 + ^ h . Analogously, Fig. 6c allows one to 
define C triples: 

r-OO />OG 

F 3] h{k m ,ti,t 2 ,t 3 ) = / dt' 2 dt' 3 T^\k m ;ti,t 2 ,t' 3 )F 2 . h (k m ;t' 2 - t 2 )F 2 . h (k m+1 ,t 3 - t 3 ) 
Jo Jo 

r-OG />OG 

+ / dt[ / dt' 3 T A . h {k m \t' 1 ,t 2 ,t' 3 )F 2 -f l {k m -\,t' 1 — t{]F 2 -h{k m+ i,t' 3 — t 3 ) 
Jo Jo 

poo poo 

+ / dt[ / <Ti+ 1 \fc ro ;4,t 2 ,t 3 )i ? 2;/ l (A; m -i,4-ii)F 2 ^(fc m ,t 2 -i 2 ) 
Jo Jo 

poo poo 

+ / dt[T B . h {k m ,t' 1 ,t 2 ,t 3 )F 2 -} l {k m -i,t' 1 — t\) + / dt' 2 T B . h (k m , t\, t 2 , t 3 )F 2;h (k m , t 2 — t 2 ) 
Jo Jo 

poo 

+ / dt l 3 T B % lh (k m MMA)F2:h{k m +iA-h)+T C ;h{k m -M 1 t 2 M) ■ (56) 
Jo 



11 



B. Representation of the 4-point functions 

This theory has four different fundamental 4-point sta- 
tistical objects: the correlation function F^h and thre e 
non-linear Green's' functions: G^^-h (jAjj), Gz^h ( A10 ), 
and G\ 3.^. These objects are shown in Fig. a as circles 
with four tails: four wavy tails for F A . h , three wavy and 
one straight tail for G^^-hi etc. As the 3-point objects, 
also the 4-point objects can be represented as combina- 
tion of lower order 2 and 3-point objects, and in addition 
they have a piece that cannot be so represented. Consider 
Fig. 7. For every 4-point statistical object we show the 
part obtained from 2-point objects, which is the Gaussian 



decomposition, a part that is given in terms of 2-point 
objects and 3-point vertices, and finally the part that 
calls for the introduction of 4-point vertices which are 
represented by empty squares and denoted by D. There 
are four different D- vertices, depending on the nature of 
the tails connecting them. As in the case of the 3-point 
vertices, the 4-point vertices that can be solved for in 
terms of the 4-point statistical objects. The figures of- 
fer a strong motivation for the graphic representation. 
Analytically the equation for F^h consist of 63 terms, 
Eqs. for G^^h has 31 terms, etc. The graphic notation 
is sufficient for all our future considerations. 



a- -V 



b. g 3 . 



d. G 2 ,, 



3 (^AAATl +45 { 



, r aaa — , r 

J I AAAW- ) +21 \ 



3 {X= } * 9 { 




FIG. 7. Compact representation of the 4th-order 
statistical objects in terms of triplices and irre- 
ducible 4th order vertices (empty squares) . Double 
lines denote either wavy or straight lines. Thus a 
short wavy and short double line denote either pair 
correlator (long wavy line, when the short double 
line denotes a wavy one) or the Green's' function 
(short wavy and short straight line, when the dou- 
ble line denotes a straight one). One understands 
a long double line as consisting of two short dou- 
ble line. Therefore a long double line may denote 3 
versions of 2-point objects: the 2-point correlation 
function and 2 orientations of the Green's functions. 
The long double line is not allowed to be made of 
two short straight lines. The different topologies 
are counted in these figures, for example, Panel a 
one has 3 different Gaussian decompositions of F^-h, 
45 different contributions in terms of triplices and 
15 contributions in terms of 4-point vertices. 



In particular, by substituting these representations on 
the RHS of the hierarchy of equation for the 3-point ob- 
jects (Panel c and bottom lines of the Panels a and b 
of Fig. H) we get Equations for the 3-point objects in 
terms of the 2, 3-point objects and 4-point vertices. An 
example of such an equation is given in Fig. 8. 



C. Closure on the level of 4th order objects 



The simplest closure which may already predict 
anomalous scaling is very simple, and is obtained by ne- 
glecting the empty squares in all the equations of motion 
for 2- and 3-point objects. Below we show that in this ap- 
proximation we account exactly for all terms up to 0(e ) 
and neglect terms of 0(e 6 ) and smaller. The resulting 
system of equations can be exactly reduced to a system 
of four evolution equations for F2-h and the three trip- 
lices. A simplified version of these equations involving 
the A-triplex only is shown in Fig. |5|. 



Equation for Q w . , 




~s/w — *- 




FIG. 8. Equation 
for G\ t 2\h- For sim- 
plicity we did not dis- 
play here the terms 
with B- and 
C-triplices. Topolog- 
ically they look like 
the exhibited terms 
with the A-triplex. 



D. Closure on the level of 5th order objects 

In order to formulate this closure consider four sets of 
hierarchical equations: 

• The first set consists of 3 evolution equations for the 
correlation functions F2-h, F^-h and F^h- These equa- 
tions contain on their RHS Fa-u, F^-h and Fs-h respec- 
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tively. 

• The second set consists of 3 equations for the Green's' 
functions Gi.i ; h, G^.i;^, and G^h, having on their RHS 
the G 2 ,i-h, G 3t i- h , and G i: i-h- 

• The third set consists of 2 equations for the Green's' 
functions G\^h and C?2,2;/i in terms of G^^h and G^^h- 

• The fourth set has just one equation for the Green's 
function G\^ t h in terms of G2,A\h- 




} = 6{ 



FIG. 9. Exact equation for Dn,&—n (empty square) 



The new objects that appear at this level are the four 
5th order functions F^-h, Gz&h, G^^-h and G^\-h which 
appear in the RHS of the last equation in each of the four 
sets. Like the 4th-order objects they may be presented as 
a sum of reducible and irreducible parts, the latter being 
related to a new set of 5th order vertices that can be rep- 
resented graphically as empty pentagons. The reducible 
part is a sum of all the topologically allowed contribu- 
tions made from lower order objects: F^h, G\^h, three 
triplices and four 4th order vertices D^ ntn -h- The ir- 
reducible part of the statistical objects is, by definition, 
what is not reducible. The Closure on level of 5th-order 
objects is obtained by disregarding all the irreducible 
contributions in 5th-order objects. 



We are not going to study explicitly the equations ob- 
tained from the closure at the level of the 5th order ob- 
jects. We use this level to justify the relative smallness of 
what was neglected in the closure on the level of the 4th 
order objects. To this aim we derive now an exact equa- 
tion for the empty squares that does not contain any time 
derivative. This is done as follows: Compute dF^h/dti 
using the representation of F^-h via Z?4_„ in; ^ and 2nd or- 
der objects, shown in Fig. 7. For the time derivatives 
of the 2nd order objects use the corresponding evolu- 
tion equations involving F^-h and the third order Green's 
functions. The latter in their turn may be expressed ex- 
actly via 3rd order vertices (black triangles). On the 
other hand we have the evolution equation for dF^/dti 
which relates it to F^-h- Now F$-h can be presented as 
a sum of reducible contributions (made from lower order 
objects) and an irreducible part, which is presented via 
empty pentagons. Equating the formulae for dF^/dti 
one gets after some tedious algebra an exact linear equa- 
tion for the 4th order vertices in terms of lower order 
objects and the 5th order empty pentagons. This equa- 
tion is displayed in Fig. |[ The closure on the 5th level 
amounts to discarding the empty pentagons, leading to 
a linear inhomogeneous equation for the empty squares. 
This equation serves for the order of magnitude estimates 
presented in the next section. 



VII. ORDERS OF MAGNITUDE 

In this Section we explain how the small parameter e 
controls the sequential steps in the closure procedure in 
the (N, e)-model. In particular, we address the empty 
square in Fig. 8 that is being discarded in the closure on 
the level of 4th order objects, showing that it is of lower 
order in e compared to the terms that are retained in the 
limit N — > oo. 

First, recall that initially all our objects had "Fourier 
space" indices a, (3, 7, . . . and spinor indices a, a 1 , a" . . .. 
The tensorial correlation functions FSt„ „, , F.f' Q and 

F2'co' are defined by Eqs. (|2^)- (|5|)and the tensorial 
Green's functions of the 2nd, 3rd and 4th order are de- 
fined by Eqs. ([A§)-(|A|). Equations @- @ introduce 
the scalar objects F2, F3 and F4, as sums over tensor in- 
dices. For these objects we derived scalar equations and 
that were represented diagrammatically. Of course, all 
the diagrammatic equations (including equations on an 
/i-slice) discussed so far may be understood also as equa- 
tions for the tensor objects with explicit a, /3, 7, . . . and 
a, a', a" . . . indices. In this case each equation (or di- 
agram) for the tensor objects includes summations over 
dummy Fourier and spinor indices for all the interven- 
ing objects, (vertices, triplices, etc.), in the body of the 
diagrams. This calls for definitions of tensorial triplices 
and fourth and fifth order vertices. In particular, each 
triplex is understood as a sum similar to (p^) of objects 
depending on a, /3's. For example: 



a, 0,1 



+ A Q , i/ 3 +7+ jv]T 4 ft a '' 3 ' 7 (fc rl ; tl, #2) £3) 



(57) 



where the sums on cr's have already been performed. The 
objects T°;f ' 7 are obtained from tensorial equations iden- 
tical to j52]-p6|), but with the 2nd and 3rd order objects 
having the appropriate Greek indices. Having done so 
is of order 0(1/ VN): this stems 



it follows that T£f ' 7 



from the fact that F^.f ' 7 is of order 0(1/ y/N) and F$ h 
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and G"°. h are O(l). In other words, T^.'f n has the same 
order magnitude and the same tensorial structure as the 
3-point correlation function F^.f' 1 . Again, Eq. ( |57| ) can 
be inverted to give: 



A-h > 



(58) 



where is a function O(l). Similar relationships may 
be written for the other triplices. 

Using this information we can discuss the order of mag- 
nitude of the terms in the equation for G\^h exhibited 
in Fig. 8. All the following considerations apply similarly 
to the equations for F^h and C?2,i;fc- In fact the order of 
magnitude of the different terms in the equations does 
depend only on the topology of the graphs and not on 
the kind of 2-point functions. 

The equation in Fig. 8 has three types of contribu- 
tions. The first one on the RHS has a "simple" loop 
(made of 2 lines representing 2nd order objects). The 
next three terms have a loop made from 3 lines ("trian- 
gular" loop) and the last term has a 4th order vertex, an 
empty square. 

The term with the simple loop has two sums over Greek 
indices, each having a product of two complex conjugate 
<&'s (three factors $'s coming from the three vertices and 
one coming from the definition of G\.2-h)- Accordingly 
this term is of the order of unity, as the LHS. In the tri- 
angular loops we also have two sums on Greek indices of 
four factors of <&. In this case the term will be multiplied 
by the following factor 



■P-6 



(59) 



a/3 8 



= ^ ^(e + Vl-e 2 ^ a+ ^){e + y/l-e 2 * a ' s > a - s ) 



af3S 



x (e + y/l-e 2 ^- S ' fi+s )(e + y/l - 6 2tfa-*.a+/S,-/J-«) _ 

The four factors \& do not appear as couples of complex 
conjugates. As a result, after averaging over the ran- 
domness and sending N — > oo all the terms that present 
a random phase in (]59| ) will go to zero and only the factor 
e 4 will survive. The conclusion, which is of fundamental 
importance for the solution of our model, is that terms 
with triangular loops are smaller than terms with a sim- 
ple loop by a factor o/0(e 4 ). This is the first place that 
e appears as an important parameter that distinguishes 
the order of magnitude of competing terms. We reiter- 
ate that usually, without the introduction of e, all the 
diagrams that one obtains in the context of turbulent 
statistics are of the same order. 

The second important role of the small parameter e is 
to make the last term with an empty square even smaller 
than the terms just analyzed. We show that it is of the 
order of e 6 . In order to see this, consider the equation for 
the empty squares exhibited in Fig. The irreducible 



part of F 5 -h is related here to the empty pentagon. At 
this point we discard the empty pentagon without much 
ado, and find the order of magnitude of the empty square. 
Of course, this step should be justified by estimating next 
the order of the empty pentagon, which can be done by 
going to the closure on the level of 6th order objects 
[with the result that the term with the empty pentagon 
is smaller than the other terms in this equations by a 
factor of 0(e 2 ). This term is of 0(e 8 )]. 

After discarding the empty pentagon we have an equa- 
tion for the empty square in terms of itself and of lower 
order objects. This equation corresponds to the top line 
in Fig. 0. On the RHS we have objects with four vertices 
in one loop. Accordingly they are of 0(e 4 ). The LHS 
has an operator that acts on our empty square. The op- 
erator contains unity, and is therefore estimated to be 
of the O(l). This is in fact a dangerous point. We are 
not guaranteed a priori that the empty square is not a 
zero mode of this operator. We are going to make an 
assumption that this is not the case. This assumption 
cannot be justified before the end of the calculation. In 
principle, one needs to proceed under this assumption, 
compute explicitly the 2nd and 3rd order objects, and 
justify the assumption a posteriori by a direct substitu- 
tion and calculation. 

Making this assumption we conclude that the term in- 
cluding the empty square in Fig. 8 is of 0(e 6 ) and it 
can be neglected with impunity. This leaves us with the 
equation represented in Fig. ^. Previously we were satis- 
fied with this graphic representation but now we need to 
discuss its analytic form. In Appendix [b| we display the 
algebraic form in the ^-r epre sentation in full. For the 
scale invariant functions ( B18 ) these equations may be 
written as follows: 

- [in + d h {n)]T 2 ;h{v) = B h (si) , (60) 
v{in + + a*(n + )}f£ h 1) (n_,n ) (ei) 

= e 4 [pi +1) (f7_,fi )+4 +1) (fi_,f7 )], 

- [ifi + a h (Sl )]7j®(Sl-,Sl + ) (62) 

= e 4 [Pi 0) (fi„,f7 + )+4 0) (^,f7 + )], 

- [in_ + a{n-)]f^ h 1 \n 0l n+) (63) 

Here Ti-^h is a dimensionless 2nd order correlation func- 
tions, "T^h are dimensionless triplices. The term with 
the mass operator a originates from the diagram A in 
the equation for Ti h, and from diagram C in the Eqs. 



for 7f) h . The terms £ h 



(0 



and T>ff\ with a factor e 4 



front, originate from the triangle diagrams D and E in 
Fig. 8. Analytical forms of these functions are given in 
Appendix [b]. 
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VIII. 



PRELIMINARY ANALYSIS OF THE 
SOLVABILITY CONDITIONS 



Equations (|60[)-(63), which are exact in the limit e — » 
0, are nonlinear integral equations for one function of one 
frequency ^2;/i(^) and three functions T^f? (for i = and 
t = ±1) of two frequencies. A complete analysis of these 
equations is very difficult and beyond of the scope of this 
paper. In this section we will only demonstrate how the 
solvability conditions of these equations may determine 
the function Z(h). For the purpose of this demonstra- 
tion it is sufficient to restrict ourselves to the 11 one-triplex 
reduction" of Eqs. (|60|)- (p3"|). Considering only the two 
equations for JF 2; /j and T\.l, we substitute on the RHS 
of these equations instead of Tjffi : 



-[»n + * fc (n)]^ a! h(n) = flh(fi), 
[m + a h (n )]f^l(n_,n + ) 
m 



(64) 
(65) 



= e 4 [^(n_,n + )+^(Q_,n + )] . 

Clearly, neglecting the difference between triplices with 
different £ is an uncontrolled step which can (and will) 
be improved in future work. 

Given a set of integral equations (|64| ) , (|65| ) it is custom- 
ary to expand the unknown functions in an appropriate 
complete set of functions, and to project the resulting 
expanded form on each function in the set separately. 
In this way one generates a countable set of algebraic 
equations. The least automatic step in this procedure is 
the choice of the complete set of functions. By choosing 
the low order functions in the basis set to represent in 
some sense the properties of the expected solutions, one 
can hope that a truncated set may serve as well. It is 
of course also very helpful if the used set of functions 
is simple enough to allow as much analytic integration 
as possible. To find a good compromise between these 
requirements we suggest to use a set of functions: 



/o(O) 



Z7 



/«(«) 



1 



g?7" Q + «7 



(66) 



which correspond in t-representation to the set of 
Laguerre functions: exp(— jt) , . . . t n exp(-jt) . . . With 
these functions one may compute analytically all the 
needed integrals and they have a reasonable asymptotic 
behavior for large t. 

The explicit solution of Eqs. ([34|), (|65|) is not a triv- 
ial task and considerable amount of work will be called 
for in the future to find and to analyze it. In this pa- 
per we proceed on qualitative grounds to demonstrate 
that the solvability condition of these equations contains 
the phenomenon of anomalous scaling. For this purposes 
we will employ just the first non-vanishing term in the 
set (pi) as displayed below: 



£a,/.(n) F 



n 2 + 7 2 ' 



(n+-i 11 ) 2 (n_ + i 11 y 



(67) 
(68) 



One sees that for T2;h(Sl) we took lm/o(0) as required 
by the evenness of a real function JF 2: ^(f2). In t- 
representation the triplex Tj^h{t—,t+) must be zero if 
t + and/or t_ are less or equal to zero. The homoge- 
neous evolution equations for the triplices dictate that 
they are small for small t+ and/or i_. Therefore there 
is no contribution of /o(^) to the triplices, and we em- 
ploy the representation ( pq) in which the choice between 
/i(f2) and /-j* (fi) is dictated by causality. We introduced 
in Eq. (JgsJ) an additional freedom (71 ^ 7) which will 
allow to see later how qualitative results depend on the 
choice of the trial functions. The results are only weakly 
sensitive to the ratio 71/7. 

As one knows for example from quantum mechani- 
cal calculations of energy levels, eigenvalues are usually 
much less sensitive to the choice of trial functions than 
the corresponding eigenfunctions. This provides an addi- 
tional justification to the expectation that our eigenvalue 
problem of finding 2(h) from the solvability condition of 
Eqs. (p^), ( pHf ) is still meaningful even with very simple 
trial functions (p7|), (pq). To formulate the solvability 
condition we consider a consequence of Eq. ( |65[ ) which 
may be called their "one-frequency" reduction: multiply 
Eq. d§) for T h (0) by S(fl + f2_ - and integrate over 
fi+ and f2_. This gives: 



[ifi + a h (Q)]f h (Q) = e 4 [f> h (Q) +4(0)] 



(69) 



in which 



f h (fl) = J da - da+s ^— a+) Ti 0) (Sl_,Q + ) , (70) 
V h (Q) = J ^- <m+J g f0 --° +) pP )) (n_,n + ) , (71) 

s h (fi) = / "+)p g(o)(n_,n + ) . (72) 

Equation ( |69| ) has exactly the same factor [iQ + ah(Q)) 
in its LHS as Eq. @ for T 2 -M n )- 



[tn + flffc(n)]^ 2 ;/i(n) = Bh(n) 



(73) 



This allows the elimination of the non-homogeneous 
terms cx Q, multiplying Eq. ( |7^ ) by T^Vt) and Eq. ( |69| ) by 
■^2;ft(0) and taking their difference: 



f h (n)B h (n) = e 4 AMn)[i > h(n) + £ h (n)} 



(74) 



Recall that according to (B24) Bu is a form linear in the 
triplices T^h and quadratic in ^2^, while T>h and T>h 
are forms quadratic in Ta\H and linear in T^n- Therefore 
Eq. ([74]) is a homogeneous form which is quadratic both 
in "Ta\H and Ti;h. With our representation (|67|), ( |6^ ) of 
T2-h and Ta-h Eq. (ff^) is proportional to F 2 T 2 which 
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may be canceled if there exists a non-trivial solution of 
the problem. Therefore one may think that Eq. ( f74| ) is 
the solvability condition which we are looking for. How- 
ever, the trial functions ( |67| ) and ( |68[ ) are not a solution 
of Eqs. ( |73|), (|69|). Therefore we cannot expect to sat- 
isfy Eq. ([74J) for all values of the frequency. According 
to our calculational scheme we have to project the ini- 
tial Eqs. (p^), ( f73| ) on the (truncated) set of functions 
fo(£l), fi(fl) and consider the algebraic consequences of 
this step. Actually Eq. ( |74| ) sheds light on how we may 
proceed in order to find a solvability condition in our 
scheme. Namely, we may project Eq. (O) for T2.i1 on 
the triplex trial function ( |68| ) and Eq. ( 63 ) for Ta,h on 
the trial function ( |67[ ) for Ti;h an d subtract the results. 
This procedure corresponds to Q- integration of Eq. (^3) : 



fe(n)B fc (n) 

Z7T 

4 ^ ^[^ 2! /,(n)© h (n) 



(75) 



•^ 2i fc(n)4(fi)] 



In Eq. ((rj) B h (fi) is given by Eqs (Egg ) - (|B27| ), V h {9) 
and fh(fi) are defined by Eqs. (|71|) and ([71|) in which 



,n+) and 4 0) ( fi - 



0+) are given by flB28| ) - 
( B30| ). All these equations contain the functions ^2;ft,(^) 
and 7/j(i7) for which we will use the representations 
@ and pi). Finally Eq. @ yields: 



i? 



-U a (h) + bU b (h) + cR z U c (h) = 0, 



(76) 



where i? 2 = A 2 +^(' 1 ) and £7 a (/i), f7 5 (ft) and [/ c (/i) are 
^-integrals of U a -h(^i), Ut-h(^i) and J7 c . h (f2) and 



I7 a;h (fi) = K;h(«) - e 4 W a;?l (n) , j7 6;h (fi) = V b , h {VL) 



U Cih (n) = V c;h (Q) - e 4 W c , h (Q) 



(77) 



Equation ( f76] ) are quadratic with respect to i? z — 
y2+z{h) an( j a ]] ows to find Z as a function of ft,. The 
actual calculation has been done with the help of Matem- 
atica, a system for doing analytical calculations by a com- 
puter JT3] . The resulting set of functions Z(K) for e 
0. 0.25, 0.5, 0.75, 1, for the spacing parameters A = 2, 4 
and the parameter b = 0, —0.5, —1 is shown in Fig. fL0[ 
Figure |ll| displays the functions Z[h) for the same set of 
e, b = -2 and A = 2, 4, 8. Note that at for e = (solid 
lines) in both Figs. 10 and 11 at h = 1 all the functions 
Z(h) are equal to —2 (this is relatively easy to see ana- 
lytically). For h < 1 one sees that the functions Z(h) go 
down almost as straight lines. An important conclusion 



is that at e = the functions Z(h) are always negative 
(moreover, Z(h) < —2). Recall that the scaling expo- 
nents Cn are related with Z(h) via saddle-point require- 
ment ([49j) and, in particular, Co = min^ Z(h). By def- 
inition, the zero-order structure function is 1, therefore 
Co = and, consequently, for physical solutions Z(h) > 0. 
Negative values of Z(h) are unphysical like the negative 
values of masses which may appear as a solution for equi- 
librium of a mechanical system. Remember that in the 
derivation of the closure equations we have neglected in- 
homogeneous terms with normal K41 scaling with respect 
of terms with anomalous scaling; this is only possible if 
Z(h) > 0. Therefore for negative values of Z(h) the 
"anomalous" solutions must be disregarded in favor of 
solutions with normal scaling. The first important con- 
clusion is therefore that at e = the qualitative analysis 
of the controlled closure predicts K^l scaling. 

We next consider the parameter range for which ac < 
0. This mimics the Navier Stokes case in directing the 
cascade from large to small scales. The parameter b is 
chosen in the region of stability — 2 < b < 0, for which 
the effect of the flux of the second ( "helicity" ) integral of 
motion becomes irrelevant deep in the inertial interval. 
With a taken to be unity, these two constraints limit b to 
the interval [—1,0]. Indeed, in most previous simulations 
the value b = —0.5 was chosen. The calculated functions 
Z(K) for this region of b are exhibited in Fig. [h]. One 
sees that in all the panels the values of Z(h) increase with 
e and some of the functions Z (h) become positive in the 
left part of the h interval. Moreover, Z(h) gains a min- 
imum (although still in the negative region). All these 
qualitative features are precisely those expected for the 
functions Z(h) in the case of anomalous scaling. We note 
that the present level of analysis cannot yield quantita- 
tive results. The functions Z[h) are still mostly negative, 
and we only see the good trend of creating a minimum 
and inching up toward positive territory with the increase 
in e. Although we cannot discuss at present the actual 
values of the anomalous scaling exponents, their depen- 
dence on the parameters of the model: e, A and b, etc, 
we think that results in the Fig. |l(] can be considered as 
definite evidence for the birth of anomalous scaling with 
the increase of the value of e. An additional encouraging 
result is the finding that in the non-physical interval of 
b (see in Fig. [ll|) we do not see any of the good phe- 
nomena discussed above: at e = we have two negative 
branches of Z[h) and with increasing of the value of e 
the solutions stay in between these branches without a 
tendency to increase. Moreover for e > 1/2 the solutions 
Z(K) disappear altogether. 
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FIG. 10. The function Z(h) for different values of e, b and A. Different line types denote different values of e: e =0, 0.25, 0.5, 
0.75 and 1 are shown with solid, long dashed, dashed, dash-dotted and dotted lines respectively. The three columns correspond 
to b = 0, —0.5 and —1, as shown over the panels. The upper panels correspond to A = 2, the lower ones to A = 4. The ratio 
7i/7=0.5 . 
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FIG. 11. The functions Z(h) at b = —2 and 71/7=0.5. The line types correspond to those in the previous figure. The values 
of A are shown over the different panels. 



IX. SUMMARY AND DISCUSSION 

For the sake of clarity we summarize verbally the im- 
portant steps in this paper. 

1. We first introduced the (N, e)-shell model in which 
N copies of the "Sabra" shell model are coupled 
with the interaction amplitudes shown in Eq. ( 2C|) . 
These amplitudes contain a deterministic and a 
random part with the parameter e determining 
their relative importance. 

2. For e = the model reduces to the well known 
random coupling model [Q, in which there is no 
anomalous scaling in the limit N — > 00. For e = 1 



the system boils down to N uncoupled copies of the 
"Sabra" model with anomalous scaling for each. 

3. We studied numerically the behavior of the model 
for values of e in the interval [0, 1] in the limit 
N — ► 00. The results are shown in Figs.l and 2. For 
e = we find no anomalous scaling as expected, for 
e = 0.8 we find anomalous corrections smaller than 
those found for the Sabra model which pertains to 
e = 1. 

4. We derived the linear hierarchy of evolution equa- 
tions for the correlation and response functions of 
this model. In the inviscid limit the hierarchy of 
equations for the correlation functions is linear and 
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homogeneous. 

5. We pointed out that the homogeneity of the hierar- 
chy of equations results in a rescaling symmetry of 
these equations which foliates the solutions on so- 
called /i-slices, where h is the rescaling exponents of 
the velocity variable. The nth order velocity corre- 
lation function rescales on an /i-slice with a scaling 
exponent nh + Z(h), where Z(h) is associated with 
the rescaling properties of the probability measure. 

6. The full solution is written as an integral over con- 
tributions with different h, and leads to multiscal- 
ing via the saddle point calculation of the corre- 
sponding integrals over h. The values of the scal- 
ing exponents are determined by the function Z(K) 
according to Eq. (ff9|) . 

7. The hierarchical equations written on an h slice can 
be closed in a controlled fashion using e as a small 
parameter. We derived the closed equations for the 
2nd and 3rd order objects, retaining terms up to 
0(e 4 ) and showing that the parts neglected are of 
0(e 6 ). We could go a further step in the closure 
scheme by writing down a system for the 2nd, 3rd 
and 4th order objects, retaining terms up to 0(e 6 ) 
and neglecting terms of 0(e 8 ) etc. 



8. Section VIII presents a preliminary analysis of Z(h) 
from solvability condition for the closed set of equa- 
tions for the 2nd and 3rd order objects. The full 
analysis is technically involved and calls for fur- 
ther work. The main qualitative conclusions of the 
present analysis is that there exists no anomalous 
solution for e = 0, whereas for larger values of e 



we see the birth of anomalous scaling. The level 
of precision is currently not sufficient for quanti- 
tative comparisons with the numerical simulations, 
but we judge the results extremely encouraging for 
further work in the near future. 

In summary, we presented above what appears to be the 
first attempt to compute from the equations of motion 
the anomalous exponents of the correlations of a nonlin- 
ear turbulent system in a controlled fashion. The small 
parameter e has been used to close the infinite hierarchy 
of equations for the statistical objects, from which one 
can compute the 2nd and 3rd correlation functions expl- 
citly. Nevertheless we considered the calculation of Z(h) 
only, and also this calculation turned out to be too hard 
to be performed exactly. The bottom line of this analysis 
is that (i) the anomaly of the scaling exponents is ana- 
lytically predictable by the method developed above, and 
(ii) the actual calculation of numbers is rather difficult 
and calls for additional careful work. 
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APPENDIX A: SKETCH OF DERIVATION OF THE HIERARCHY OF EVOLUTION EQUATIONS 

1. Hierarchy for the correlation functions 



Let us compute dF%. a a {k n \t — t')/dt defined by Eq. ( |23| ) at t' — 0, substituting instead of du" a (t)/dt the RHS of 
Eq. (|l9|). With the help of definition (|2~i| ) the result may be written as: 



dt 



AZ,„„ 



x {ak n+1 F™f;? a „(k n+1 ;0,t,t) + bk n F«f^„(k n ;t,0,t) + ch^F^^n-M,®)} 



(Al) 



where we used the symmetry properties A" a i a „ = A^.^,, . Summing this equation over a and er we derive Eqs. (pq ) with 
the invariant combination (B4h in the RHS. 



2. Definitions of the Green's function 



aa:defG 



First we define the "tensorial" Green's function 
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(A2) 



and its scalar counterpart which is invariant under all the symmetry transformations described in Subsect. IV A: 

G 1 ^k n -M-t 2 ) = ^Y. G ^A k -M-t2) ■ (A3) 



Next we define the 3rd order tensorial Green's functions G2.1 and G1.2 describing a response of two velocities on one 
forcing and vice versa: 



rW3,7 ft, u Lit + \- I J ["m. g '( f m)M?; g »fa)] \ 



G\ 2:aa' a" i^n: ^"mi kt\ En, ^m: ^l) 



(A4) 
(A5) 



Analogously to the definition of the scalar 3rd order correlation function ([34|) we define here the scalar counterparts 
of these 3rd order Green's functions: 

a, 



a,/3,7 a : a' a" 



^ ^2 l-cra' a" i^ni &mi ^nj ^W) ^) ; 

G\^{k n ,km, ke; t n ; t m: ig) = — ^ ^ A a a , a ,> [A a ^ +7 + A a+ jv : /3+7 + A a;/ 3 +7+A r] 



a,/3,7 a,o' a" 
-ra,/3,7 



(A6) 



(A7) 



As we discussed in Subsect. VIA any three point objects of the Sabra-shcll model differ fro m zero only for three 
consecutive shell numbers. This allows one to consider instead of each Green's function (A7) and (A6) of three k 
arguments three functions with just one (middle) k argument. Exactly like in the case of triplices we introduce 
superscripts and ^ to distinguish them. Using the convention discussed in Subsect. VIA we can define G^l, 
, G^l and G21 • For example: 



G2 1 {km ? im—1 5 Em j ) — G2.I (km 5 km— 1 5 5 tm 5 tm— 1 ) ^m+l) ■ 

Similarly we define the successive scalar response functions of the 4th order as follows: 

<J3,l(#nl, K n2, Kn3, /s n4; El, S2j E3, ti) = J^p A a " a 3 A a t a & \ 



1 



^(^, CT »fe)<3, CT3 («3))' 



G 2 ,2(fc™i, k n2 ,k n?n k n ^ ti, t 2 , t 3 , U) = -tt^AI„ a3 A° ias . — — g- 

JV \ <5/„l,<T'( I l) (5 /n4, ( T 4 (E4) 



(A8) 



(A9) 



(A10) 



with the convention of sums over repeated indices. Note that all the G's are real functions. Again we note that it is 
more convenient to write the hierarchical equations in the terms of quantities summed over cr's. These are invariant 



under all the symmetry transformations discussed in Sect. IV A 



3. Hierarchical equations for the Green's functions 

The tensorial equation for the 2nd order G"'" (T(T reads: 



dt 



G l';i-a,A k ^ f ) = E $7 ^'"[ A «+/3,7 + A Q+/3+ Ar, 7 + A a+fj ^ +N ]{A°, a „ [la.n+lG^ al<jll {k n , k n+ l, k n+2 ] 0, t, t) 



{3,7 , a' .a" 



+7&,»G^;7 CT , CT ,,(fc n ,^ . (All) 
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By following more or less the same procedure as for the correlation functions in Appendix A 1 we get the hierarchical 
equations: 

— Gi,i{k n -t)= la , n+ iG^\k n+1 -0,t,t^ . (A12) 



Analogously for the higher order Green's functions one has: 

-G 2 i(k n ; tl, *2, £3) = 7a,nG ! 3 i i(fc n , k n , fe n +ife n +i; f2, t\,tl,tz) + 7fe, n -lG ! 3 ) l(fc J j, fc n _2, fcn&n+ljfe, tl, tl, ts) 



d 



+7c,n-2G , 3,i(/c„, fc n -3, fcn-2fcn+i;*2| t\,t\,h) , (A13) 
-G^ 2 ■* (^ni ^1) *3) = la,nG-2,2 (k n , k n +l, k n , k n+ i, tl, tl, t%, t$) + 76^-1^2,2 (^n-2i ^n, k n k n+ i;ti, tl, t 2 , £3) 

+Jc,n-2G 2 ,2(k n - 3l k n - 2 , k n , k n+ i;ti,ti,t2,t 3 ) , (A14) 
and similarly for the other G2 i's and Gi 2's. 
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APPENDIX B: THE CLOSURE EQUATIONS IN 
THE A-TRIPLEX APPROXIMATION 



1. Equations in the ^-representation 



Given a stationary case, when all the functions depend 
on time differences only, it is useful to consider the clo- 
sure equations in the cj-representation. Define: 



F 2 -h{k n \ w) = / F2- h (k n ;t)exp(iujt)dt, (Bl) 

— OO 
OO 

F* M {k n ;uj)= [ F* h (k n ;t)exp{-iujt)dt . (B2) 



As in these equations, we will use tilde signs over a char- 
acter to denote a Fourier transform of that function. It 



follows from Eq. (Bl) that 



F 2 (k n ;t) = / F 2 (k n ; lu) exp(-iu;t) 



2^ 



(B3) 



According to the definition ( |2^j3^ ) F 2 -h(k n ;t) — 
F^-hikn; — t). Therefore F2 ; h(k n ;uj) is a reai function of 
to. 

Next we define the Fourier transform Th of any 3rd 
order object Th as follows: 



27r7)j(fc„; luq, lu + )5(lu^ + luq — lu + ) = j dt-dtodt + 

—00 

xT h (k n ;t-,t ,t + )ex.p[i(Lj-t- +aj t Q - uj + t + )} . (B4) 



Recall that with our convention (see Subsect. VIA) to 
and t± are times of the n- and (n ± l)-shells. With this 
definition the sum of incoming frequencies cj_ + luq is 
equal to the outgoing frequency uj + . 

The equations for the 2nd order correlation function 
(two top lines in Fig. ||) and for triplices (two last lines) 
in the tj-rcpresentation may be written as: 



- [iuj + Y, h (k n ;u)]F2-h{k„;uj) = B h (k n ,u>) , 
[iw+ + S^(fc n ;w + )]f)[ +1) (fc n _i; w_, w ) = e 4 [D h (kn-V, w ) + E h +1 \k n -i; w )] 

— [iu + Eh(k n ]Wo)]f^(k n ',U-,tjj + ) = e i [E ( h 0) (k n ,uj^,uj + ) + D h °\k n ,uj-,LU + )] , 
[iu- + S/ l (fc„;w_)]T , ^ 1) (fc n+ i; u , uj + ) = e 4 [D h ~ 1 \k n+ i;coo, w+) + E h ^ 1 \k n+1 ; u)q, 



(B5) 
(B6) 
(B7) 
(B8) 



Note that the parameter e appears explicitly here, as it was discussed in Sect. VII. Here the "mass-operator" E/ t (fc n ; u>) 
originates from the diagram A in Fig. ^ and may be written as: 



T,h(k n ;u>) = — / duj-dujQduj+5(uj- + luq — lu+) 
2n 



(B9) 



Xi7n+1 -1 



F 2; h(k n+ 2]u+)T^ h *(k n+ i-,LJ-,u> + ) + F 2 - t h{k n+ i;uJo)T£ h 1 ' r (k n+ i;uj-,uj ) 5(oj - U-) 



(+1)* 
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-7n,0 



F 2 ;h(k n +i; ^+)T { A .h*{k n - 1 LJ ,uj + ) + F 2;/l (fc„_i; w_)T^ 1J *(fc„; w ) 5{u - uj q ) 



+7n-l,X 



F2;h(^„-i;w )T^. /( 1) (fc n _ 1 ;a;o,w + ) + F 2 ;/j(fc n _ 2 ; w_)T^(fc n _!;u;_, w + ) 8{u - u + )\ . 



Here and below we omit, for the sake of brevity, the frequency argument of the leg with the Green's function in A- 
tripliccs and the corresponding argument in D and E functions. It may be found from the relationship lo + = w_ +loq. 
The term Bh(k n ,uj) in the RHS of Eq. ( p35| ) corresponds to the diagram B in Fig. || 



B(k n \ui) 



2tt 



<fw_dw dw+ jln+i -iF2-h{k n+ i;Lo )F 2 -h{k n+ 2; uj + )f[.^* (k n+1 ;Lu ,ui + )6(uj - w_) (BIO) 



Functions -D correspond to the diagram D: 



i){j +1) (fc„_i; w-, w ) 



£)[ 0) (A; n ;cj_,cj + ) 



£>i 1} ( fc n+i; w o,^+) 



^{7«,or^(fc n+ i;w,w + )f^" h 1) *(fc„;w,w )A;h(fcn;w) 



+ 



"7n+l -l^AA 



* (fcri+2 ; , Lu)f^ h (k n+2 ;UJ0,u)F2;h (kn+3 ] w) j . 

Functions i? correspond to the diagrams Ea and E/3 in Fig. § 



^h +1) ( fc n-i; UJ ) 



duj f 



<j7n,o F 2;h(fcr 1 +i;w)T , j°) l (fc n _i;w_,w - uj )T^.' h (k n ; u , u) 



(°), 



2tt I 



+7n-14^2;h(fcn-2; (jj)fffi (k n -2i ^0 ~ W_ )f^* (fc„_ 2 J W, W ) 
+7n,Q F 2;h(kn-l 1 Uj)f^\k n ; UJ ,UJ + £J_ (fc n _ 1 1 U_ , w) 

+7„-l,l#2;fe(^n-i;w)^l°lr(^n-2;W - W_, Wo)Tjj^(fc n _ 2 ; W_, 



\jn+i,-iF 2 ;h(k n+2 ;uj)f'^ l l\k n ;uj-, u - u + )T [ A .^' r (k n+1 ;uj + ,uj) 



2vr L 

+7n-i4^2;/ l (fcn-i; ^)r^ l 1) (fc„_ 1 ; w + - u)f^ h (k n ; UJ,U + ) 
+ r rn+i-iF 2 -h(k n +i;uj)f^*(k n+ i;uj - u}-,u + )f^ h (kn;cj-,w) 

+7n-l,1^2;fc(fen-2;w)Tj^ 1) ^ , 

OO 

^{7n,o^2;h(fc«+i;w)f^7 ft 1) *(fc, 1 ; W + - C^WojT^^fcn+i; W,U)+) 



(Bll) 



(B12) 



(B13) 



(B14) 



(B15) 



(B16) 



(B17) 



+7 n +i,-li ? 2; h(fc„+i; w)T^) l (fc„ +2 ; oj , to + LO+)T^' r (k n+2 ;uj, u> + ) 
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2. Closure equations for Scale-invariant Functions 



All the objects appearing in the previous equations (B5)-(B17) are scale invariant, as they are defined in an "/i-slice" . 
We make use of this invariance to simplify the closed set of equations by introducing the following scale invariant 
dimensionless representation: 



fc ^0J n . h / \Un;h- / 

jj2 X n{l-h)-{h+Z(h)+X)^J_^_ \ ^ 



F 2 - h {k n \uj 
B h (k n ;uj 
T ( ^ h (k n ;Lu a ,U! b 

Here the characteristic frequency of n-shell on an /i-slice is 

uj n , h = Uk m n , /x = A 1 "' 1 









UJn;h ' 



















1= -1,0,1 



(B18) 
(B19) 
(B20) 
(B21) 
(B22) 

(B23) 



A is a spacing parameter in the shell models determined as the ratio between two consecutive shell momenta: A = 
k n +i/k n . Recall that the dependence on a single k n is represents in the case of 3rd order functions a dependence on 
three consecutive shell momenta k n _i, k n and fc n+1 . The superscript ^ determines the fc-argument of the "special 



leg" accordingly to the convention of Sect. VIA 



In the dimensionless form Eqs. (B9) are scale invariant, independent of the shell number. In these Eqs. 
B h (Q) = -?rV a]h (n) + bV b . !h (Sl)+cR z V c . lh (n), where R z = \ 2+z(h) , 



v a . h (n) 

V b;h (fl) 

v c , h {n) 



2^ ^ ' 

2^ 



fir 



dOxdO^Cg - jV- Hg) ^ m ^ , n 2yr(+i) m n ^ 



The RHS of Eq. (|62|) for T^) h has the form 



4 0) (^-,^+)+©r( fi -»°+) = ■^-W . h (fi_,n + )+cfl Jt W c . h (n_,n + ), where 

i^ k( si_.si + ) = / |J{,?# 2J ,(n)fi«>*(2,2±)fji 1 '(n_,si) 
Wc. h (n_,n + ) = | ^{i^ 2;fc (n)7^- 1) (o M ,n +/ x)7^- 1) (o,o + ) 

i^.^tj-^Q-, (r>+ - n)/i]fj°>(n,n+) + i^ a . & (n/i 2 )ri+ 1) (n / i ) n_/i)r^ 1) (n + n_,n+)} . 



(B24) 
(B25) 
(B26) 
(B27) 

(B28) 
(B29) 

(B30) 



/' 
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